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Deuteron and triton magnetic moments from NMR spectra of the hydrogen molecule 
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We present a theory and calculations of the nuclear magnetic shielding with finite nuclear mass 
effects and determine magnetic moments of deuteron and triton using the known NMR spectra of 
HD and HT molecules. The results fid = 0.857438 234 6(53) fiN and fit = 2.978 962 471(10) fiN are 
more accurate and in a good agreement with the currently accepted values. 
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When a molecule is placed in a homogeneous magnetic 
field B, its nuclei experience the field that is shielded 
by the surrounding electrons (1 — a) B. The magnitude 
of the shielding factor a, typically of the order 10“^, 
depends on the relative distance between nuclei, which 
means that nuclear magnetic resonance (NMR) can be 
used as a tool for obtaining information on the structure 
of complicated molecules. However, when the molecular 
structure can be calculated, high-precision NMR spec¬ 
tra can be used to determine the relative magnitude of 
nuclear magnetic moments [1]. Namely, the ratio of nu¬ 
clear magnetic moments is proportional to the ratio of 
measured frequencies that flip the nuclear spins 

fiA (1 - (Ja) _ /a Ia 
flB (1 — (Jb ) Jb Ib 

The accuracy to which fiB known can be transferred 
to fiA provided that one knows with sufficient precision 
the (7a — O’B difference. Here, we report on results of 
the calculation of this shielding difference for HD and 
HT molecules. Employing these results we determine 
the magnetic moment of deuteron and triton using re¬ 
cent high accuracy measurements of the proton magnetic 
moment [5] and the ratio of spin flip frequencies nn. 

Magnetic shielding of the nuclear magnetic moment 
due to the surrounding electrons has been first consid¬ 
ered by Ramsey in [7] with the help of the nonrelativistic 
Hamiltonian in the external magnetic field. His result for 
the isotropic shielding factor is 

(2) 

b 

+(<()iy]xa , 

a b b ^ 

where ex is the fine structure constant, Xb is the posi¬ 
tion of the electron b with respect to the nucleus, pb is 
its momentum, H is the molecular Hamiltonian in the 
Born-Oppenheimer (BO) approximation, and (/> is the 
electronic BO wave function. An immediate conclusion 
that can be drawn from this formula is that the shielding 
of proton and deuteron (triton) in the HD (HT) molecule 
is the same. Upon averaging with the nuclear function y. 



0-(o) = (x|cri°)(i?)|y), the difference in the shielding still 
vanishes, not only for the ground state but also for any 
excited state. Clearly one has to go beyond the BO ap¬ 
proximation and also include finite nuclear mass effects 
in the coupling to the external magnetic field. There 
have been several attempts to calculate the shielding dif¬ 
ference in HD (HT), and they are all incorrect or in¬ 
complete. Neronov and Barzakh in 1977 [1] derived the 
formula and obtained the result of 


(5cr(HD) = crd(HD) - crp(HD) = 15.0 • 10"^ 
(5cr(HT) = (Tt(HT) - cTp(HT) = 20.4-10"^ (3) 


but they started with the incomplete Hamiltonian, i.e. 
their formula (4) does not include the nuclear spin-orbit 
interaction (see 5^ — 1 terms in Eq. J0 below). Later, 
the calculations by Jaszuhski et al. [8] simulated nona- 
diabatic effects by an artificial charge difference. Their 
result of <5 (t(HD) = 9 • 10“®, although of the correct mag¬ 
nitude, is not well substantiated from the physical point 
of view, nor is it complete. Finally, in the most recent 
calculations, Golubev and Shchepkin [3] used more real¬ 
istic treatment of nonadiabatic effects, but their result of 
(5cr(HD) = 9 • 10“® was also incomplete. Undoubtedly, 
the coupling of electron motion to the nuclear motion is 
partially responsible for the difference So, but this is not 
the whole effect. 

A derivation of the finite nuclear mass correction to the 
shielding closely follows that of Ref. 13 HU] and starts 
with the Hamiltonian for electrons and nuclei, which 
includes coupling to the external electromagnetic field 
and all possible nucleus A spin-orbit interactions, i.e. 
(7i = c = 1) 
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where tt = p — eA, A is an external magnetic vector po¬ 
tential, gA is a p-factor of the nucleus A, which is related 
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to the magnetic moment by = &A dA Ia /and 
is a Coulomb interaction between electrons and nuclei. 
In order to derive a formula for the shielding constant, 
including the finite nuclear mass corrections, we perform 
two unitary transformations tp, 

H = e-^^ + dtip. (5) 

The first transformation places the gauge origin at the 
moving nucleus A. We assume that the molecule is neu¬ 
tral and that the magnetic field is homogeneous and there 
is no electric field, so 


ishes, thus pA = —PB — the independent posi¬ 

tion variables are Xa and Xb- Consider now the electronic 
Schrodinger equation 

+ = (“) 

^ a ■' 

and the next transformation ip 2 of the form 

P2 = -—^XaPB, (11) 

m. A 


= Y.^{x\A^ ^\x\xiA\^ 

a 

+e_B A*-I-- a;3 ^, (6) 


where A = A(ryi), and Xa = Xa — ta- The transformed 
momenta are 


which simplifies the nuclear kinetic energy, see Eq. (18) 
below. The electron momenta are changed to 


Pa = e =p’a 
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ruA 


Pb- 


( 12 ) 


and the potential 



= Pa + y (^a X By , 

(7) 


= Pb + ^^{^bx By , 

(8) 


1- 



= p^^ + -{Dx By , 

(9) 


where D = s Xa + eB xb is the electric dipole moment 
operator. We also assume that the total momentum van- 


V'-£ = 

« +— VfaVB(E-£), (13) 

rriA ^ 

where we omitted higher order terms in the electron nu¬ 
cleus mass ratio. The new Hamiltonian iJ after both 
transformations with pei = J2aPa- 
comes 
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From this Hamiltonian one derives the shielding defined 
by 


H,s = -7r^gAlA{l-v{R))B. (15) 
2 mA 


with 


1111 1 
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and where 


Let us consider firstly the nonrelativistic Hamiltonian 
iJnrei = H + H^, where 


Hn 



(SgPa)^ ^ m 
2 mA VIA 


yXa-VB{y-£), (18) 
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with TTin = mAViB/(viA + vib) being the nuclear re¬ 
duced mass. We assume the Born-Oppenheimer approx- 
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imation and include perturbatively using nonadia- 
batic perturbation theory (NAPT) pJl . In the zeroth 
order, the shielding is given by Eq. Q. The leading 
nonadiabatic correction to the shielding is linear in 
the electron-nuclear mass ratio, and we split it into four 


parts m- = (Tn + CTd -b tTs + O'!. The explicit for¬ 
mulae for the isotropic shielding in H 2 and isotopomers: 
CA = gb = —e, with R = tab = ta — ^b = —Xb are 
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( 20 ) 

, ( 21 ) 
( 22 ) 


where P = —i’Vfj., and {(j)\ Ar lY) = ~r 4’\'^Rp)■ 
The result in Ref. uni contains a mistake in Ud, which 
we correct here. 

In the numerical evaluation of Eqs. ( T9p2 ), the ground 
state wave function is represented with the explicitly cor¬ 
related Gaussian functions of the form 


_ g O'IAPia ^2AP2A ^2B “12 ^12 ^ (23) 


The resolvent l/(£ — PI) includes the sum of E and 11 
states, which are represented as 

4e- = P ■ {Pa X r2A) Yt.+ (24) 

(4n = P X (t)^+ . (25) 


All of the above matrix elements are in electronic vari¬ 
ables, so the derivatives with respect to the internuclear 
distance have to be obtained in advance. For this purpose 
we use the fact that total angular momentum vanishes on 
(j), so 

PxVB|(/)) = -(^a;i,xVi,) !</)), ( 26 ) 

b 


and from the P-derivative of the Schrodinger equation in 
the BO approximation we get 

= (27) 

Since we calculate only the shielding difference, all of the 
terms with the reduced mass, including the most diffi¬ 
cult Ar, are omitted. Calculations are performed us¬ 
ing 128 and 256 basis function for each symmetry with 


global optimization of all nonlinear parameters. Numer¬ 
ical results in the range R G (0, 6) a.u. are presented in 
Table |Tj Although the results for P > 3 are not in prin¬ 
ciple needed, we use them for testing against the known 
separated atoms limit, which is 

On the other hand, at small P we observe 1/P behavior, 
which is not an artifact but a result of the shielding of 
the nucleus A by the nucleus B, see Fig. 1. Our results 
at P = 1.4 a.u. in comparison to the known previous 
calculations are presented in Table |TT] 

The total isotropic magnetic shielding a is obtained 
by averaging with the nuclear wave function x, tr = 
(x|cr(°4P) -b cr*'34P)|x)! where x is a solution of the nu¬ 
clear radial equation with the BO potential augmented 
by adiabatic correction. Since the measurement [3] was 
performed ai T — 300 K we include contributions from 
the excited rotational states up to J = 9 according to 
Boltzmann distribution. The averaged result for the 
(Ramsey) shielding factor of the H 2 and isotopomers is 
= 26.33517(20) • 10“®. However, does not in¬ 
clude relativistic corrections that are of the relative order 
of ~ 10“"^. The related estimate from Ref. [5] is much 
smaller than the relativistic correction to the shielding of 
the individual hydrogen atoms, so it is possibly incorrect. 
The shielding differences in HD and HT are presented in 
Table mu Their uncertainties come from the unknown 
higher order nonadiabatic corrections. From these dif- 
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TABLE I: Shielding as a function of the internuclear dis¬ 
tance R. Results to be multiplied by 10“®. The numerical 
uncertainty of is (20), while of Sa is (2) on last digits. 


i?[a.u.] 


(5cr(HD) 

(5(t(HT) 

0.00 

59.936 77 

OO 

OO 

0.10 

58.308 74 

0.053 997 

0.069 030 

0.20 

54.933 43 

0.030 630 

0.038311 

0.40 

47.435 08 

0.020 577 

0.025 604 

0.60 

41.035 38 

0.018685 

0.023 329 

0.80 

36.014 57 

0.018680 

0.023 243 

1.00 

32.118 10 

0.019197 

0.023 663 

1.10 

30.504 37 

0.019486 

0.023 874 

1.20 

29.074 46 

0.019 759 

0.024 044 

1.30 

27.804 09 

0.019991 

0.024 153 

1.40 

26.672 56 

0.020181 

0.024 201 

1.50 

25.662 62 

0.020 324 

0.024 186 

1.60 

24.759 25 

0.020416 

0.024 107 

1.70 

23.949 94 

0.020 464 

0.023 973 

1.80 

23.224 07 

0.020 469 

0.023 792 

1.90 

22.572 23 

0.020 436 

0.023 569 

2.00 

21.986 84 

0.020 367 

0.023 308 

2.10 

21.46105 

0.020 261 

0.023015 

2.20 

20.988 70 

0.020131 

0.022 700 

2.30 

20.56501 

0.019 973 

0.022 362 

2.40 

20.185 30 

0.019 792 

0.022 010 

2.50 

19.845 67 

0.019 587 

0.021643 

2.60 

19.542 83 

0.019361 

0.021265 

2.70 

19.273 36 

0.019117 

0.020 880 

2.80 

19.034 52 

0.018855 

0.020488 

2.90 

18.823 95 

0.018 570 

0.020 085 

3.00 

18.638 83 

0.018 280 

0.019688 

3.20 

18.33714 

0.017652 

0.018886 

3.40 

18.11334 

0.016989 

0.018092 

3.60 

17.952 68 

0.016308 

0.017322 

3.80 

17.842 09 

0.015 640 

0.016598 

4.00 

17.769 75 

0.015002 

0.015 930 

4.20 

17.725 68 

0.014405 

0.015 320 

4.40 

17.70148 

0.013866 

0.014 780 

4.60 

17.690 67 

0.013 394 

0.014311 

4.80 

17.688 47 

0.012 976 

0.013899 

5.00 

17.69128 

0.012 632 

0.013557 

5.20 

17.696 74 

0.012 344 

0.013 270 

5.40 

17.703 34 

0.012 099 

0.013 023 

5.60 

17.71010 

0.011889 

0.012810 

5.80 

17.71648 

0.011 724 

0.012 638 

6.00 

17.722 24 

0.011 584 

0.012 489 

oo 

17.750 45 

0.010 753 

0.011326 


TABLE II: Extrapolation to a complete basis and comparison 
with the previous calculations of the isotropic nonrelativistic 


BO shielding in HD (HT) 
multiplied by 10“®. 

at i? = 1.4 a.u. 

Results to be 

size/ref. 


5cr(HD) 

i5cr(HT) 

128 

26.673 193 

0.020188 2 

0.0242145 

256 

26.672 556 

0.020 181 2 

0.024 2010 

512 

26.672 422 

0.020179 9 

0.024199 6 

oo 

26.672 387(35) 

0.020179 6(3) 

0.024 199 4(2) 

[T^ (1995) 

26.8139 



[T31 (1996) 

26.680 



ID (2011) 

26.677 111 




ferences and from Eq. ([^ we obtain the magnetic mo¬ 
ments of deuteron and triton, see Table |III| The ratio 
of shielded magnetic moments for HD is the average of 
two independent and consistent measurements nig and 
[3]. The obtained result for is more accurate and in 
good agreement with the presently accepted value, which 
was obtained from an unpublished experimental result by 
Philips et al. (1984), see [T3] for details. Regarding /it, 


FIG. 1: The difference 5(7(HD,i?) in ppm of the shielding 
constant between the deuteron and the proton in HD as a 
function of the internuclear distance R. The horizontal line 
is a separated atoms limit, while the dotted line is the 1/R 
asymptotics that comes from the direct interaction between 
nuclei. 

<5o-(R) 

0.06 

0.05 

0.04 

0.03 

0.02 

0.01 



TABLE III: Determination of the deuteron and the triton 


magnetic moments, /ip(HD) denotes the shielded magnetic 
moment of the proton in HD molecule. 



value 

Ref. 

Mp 

2.792 847 350(9) ^lN 

|2] 

i5ct(HD,T = 300K) 

0.020 20(2) • 10“® 

this work 

Mp(HD)/A.a(HD) 

3.257199 514(21) 
3.257199 531(29) 
3.257199 520(17) 

m 

m 

averaged 

^id = Aid(HD)//jp(HD) (1 -1- Str) fip 

0.857 438 234 6(53) mn 
0.857438 230 8(72) Atw 

this work 

nn 

Sa{HT,T = 300K) 

0.024 14(2) • 10“® 

this work 

Mt(HT)//.p(HT) 

1.066 639 893 3(7) 

1 

= Mt(HT)//.ip(HT)(l-|-5a-) 

2.978 962471(10) MN 
2.978 962 448(38) /ijv 

this work 

nn 


the CODATA [2] value is based on the earlier, less ac¬ 
curate work [5], while we use a more recent one and 
correct the value for 6a, which leads to even smaller un¬ 
certainty of /i(. We should note however, that the result 
of [B] needs confirmation and because of the lack of in¬ 
formation on temperature, we assumed T = 300 K. 

In summary, we have determined improved values for 
the deuteron and triton magnetic moments. This demon¬ 
strates that NMR spectroscopy combined with precise 
calculations of the shielding factor may lead to the most 
accurate determination of nuclear magnetic moments. 
Moreover, one may consider improving the determina¬ 
tion of magnetic moments of other light nuclei, such as 
^He, since the shielding factor for molecular hydrogen 
can be calculated including nonadiabatic and relativistic 
effects as has been already done for ^He atom Ha- 
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